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Experimental results for covalent glasses have highlighted the existence of a new self-organized 
phase due to the tendency of glass networks to minimize internal stress. Recently, we have shown 
that an equilibrated self-organized two-dimensional lattice-based model also possesses an interme- 
diate phase in which a percolating rigid cluster exists with a probability between zero and one, 
depending on the average coordination of the network. In this paper, we study the properties of this 
intermediate phase in more detail. We find that microscopic perturbations, such as the addition or 
removal of a single bond, can affect the rigidity of macroscopic regions of the network, in particular, 
creating or destroying percolation. This, together with a power-law distribution of rigid cluster sizes, 
suggests that the system is maintained in a critical state on the rigid/floppy boundary throughout 
the intermediate phase, a behavior similar to self-organized criticality, but, remarkably, in a thermo- 
dynamically equilibrated state. The distinction between percolating and non-percolating networks 
appears physically meaningless, even though the percolating cluster, when it exists, takes up a finite 
fraction of the network. We point out both similarities and differences between the intermediate 
phase and the critical point of ordinary percolation models without self-organization. Our results 
are consistent with an interpretation of recent experiments on the pressure dependence of Raman 
frequencies in chalcogenide glasses in terms of network homogeneity. 

PACS numbers: 05.65.+b, 61.43.Fs, 61.43.Bn, 64.70.Pf 



I. INTRODUCTION 

The concept of rigidity percolation, first in- 
troduced about 25 years ago by Thorpe Q 
based on work by Phillips @, describes how 
an elastic network goes from floppy to rigid 
as constraints are added to it. This theory 
has been applied with success to many sys- 
tems, including covalent glasses 0, El and 
proteins @- In the last decade, however, ex- 
perimental studies have shown that the rigidity 
phase diagram could be more complex than ini- 
tially thought, uncovering the presence of an 
intermediate phase between the floppy phase 
and the stressed-rigid phase, with the system 
in the intermediate phase bein g rigid but un- 
stressedjljl H ©, 03 El El Mill 13 13 

MM\mm 

A basic explanation for this new phase was 
first proposed by Thorpe et al. [3j ■ It was shown 
that when the network self-organizes in order 
to minimize the stress, the rigid but unstressed 
intermediate phase can indeed arise. 

In the original work by Thorpe and collabo- 
rators 0, EI , as well as subsequent simplified 
models of chalcogenide glasses by Micoulaut 
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and Phillips 22, 23], networks were constructed 
using an "aggregation" process, in which bonds 
or simple network units were added to the net- 
work without subsequent equilibration. More 
recently, Barre et al. [24| have considered a 
thermodynamically proper model with an en- 
ergy cost associated with stress and showed 
that in the canonical ensemble, the intermedi- 
ate phase still exists. In a recent paper, we have 
confirmed this result for the T = version of 
the model of Barre et al., but on a more real- 
istic regular lattice, and also have shown that 
the intermediate phase is entropically feasible 
in actual physical systems. In both the model 
by Barre et al. and our model, the intermedi- 
ate phase has an unusual property: both perco- 
lating and non-percolating networks coexist in 
the ensemble at all mean coordination numbers 
within the intermediate phase. 

While providing a general support for self- 
organization, these previous studies did not 
look in detail at the properties of networks 
in the intermediate phase. Here we provide 
a first glimpse at some of these properties. 
First, we show that in both percolating and 
non-percolating networks, the sizes of non- 
percolating clusters have a power-law distribu- 
tion. In effect, the system remains in a critical 
state over an extended range of mean coordina- 
tions, corresponding to a self-organized critical 
phase [25], but in an equilibrium system. Sec- 
ond, we find that adding or removing a single 
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bond affects the rigidity of macroscopic parts 
of the network in the intermediate phase and 
in particular, can turn a non-percolating net- 
work into a percolating one and vice versa. This 
property supports the interpretation by Wang 
and co-workers [Tsj of the puzzling response of 
vibrational frequencies to applied pressure that 
was observed in their experiments. Using our 
results, we address some intriguing questions 
that have to do with the unusual coexistence 
of percolating and non-percolating networks in 
the intermediate phase. In particular, we show 
that percolating and non-percolating networks 
can be considered identical in the intermediate 
phase, as both stay on the edge of percolation. 

This paper is structured as follows. In the 
next section, we briefly review the intermedi- 
ate phase. We then present our methodology. 
In the fourth section, we present our results on 
the properties of rigid clusters (both percolat- 
ing and non-percolating) in our model. In Sec- 
tion V, we look at the response of the network 
to local perturbations. In Section VI, we dis- 
cuss how our results can help understand the 
experiments by Wang et al. [Tg| . We finish with 
conclusions. 



II. THE INTERMEDIATE PHASE 

Using the mean-field approach first intro- 
duced by Maxwell (2(| and known as Maxwell 
counting, we can define the rigidity of a network 
in terms of its number of zero-frequency mo- 
tions, or floppy modes, F. In a d-dimensional 
network, each atom has d degrees of freedom. 
In an unconstrained network of N atoms, each 
degree of freedom corresponds to a floppy mode 
and thus F = dN. Assuming that each added 
constraint takes away one floppy mode, we can 
write 

F = dN~ N c , (1) 

where N c is the number of constraints in the 
network. As N c increases, F decreases, until 
F = is reached, and then there are no floppy 
modes left and the network undergoes a rigid- 
ity transition from floppy to rigid where a per- 
colating rigid cluster emerges in the network 
(rigid clusters are sets of mutually rigid atoms) 
Disordered networks can be conveniently char- 
acterized by the mean coordination (r), which 
is the average number of bonds connecting a 
site. In the mean-field approximation, the lo- 
cation of the rigidity transition only depends 
on (r) and not on other details, such as frac- 
tions of sites with a particular coordination. 



The transition occurs at (r) = 4 for the tri- 
angular lattice of elastic springs and (r) = 2.4 
in chalcogenide glasses; in the latter case, it is 
assumed that both bond-stretching and bond- 
bending constraints are taken into account. 

To go beyond the mean-field theory, correc- 
tions must be made. For example, adding a 
constraint to an already rigid region does not 
remove a floppy mode. Such type of constraint 
is called redundant. Redundant constraints in- 
troduce stress into the network. Such con- 
straints do not change the number of floppy 
modes (and thus violate the assumption of the 
Maxwell counting); taking this into account, 
Eq. Q becomes 

F = dN — (N c — N r ), (2) 

where N r is the number of redundant con- 
straints. Another type of correction is due to 
the fact that even above the rigidity transition 
there can still be some floppy inclusions and 
thus F > at the transition. 

To find N r , Jacobs and co-workers have in- 
troduced a topological algorithm, the pebble 
game [27l |28| . This algorithm is based on a 
theorem by Laman 29] which states that in two 
dimensions a generic network with N sites and 
B bonds does not have a redundant bond if 
and only if no subset of the network containing 
n sites and b bonds violates b < 2n — 3. A sim- 
ilar approach works in 3D, but in general only 
for bond-bending networks like those used to 
model chalcogenide glasses. 

The pebble game, described in more detail 
in Ref. [2^, characterizes the global rigidity 
of a network, provides its complete decompo- 
sition into rigid clusters and finds stressed re- 
gions. The approach uses only the topology of 
the network and not its exact geometry. Us- 
ing the pebble game, it was possible to show 
that the rigidity transition occurs at (r) = 3.961 
±0.001 [27j for the central-force triangular lat- 
tice and at (r) m 2.385 3] for an amorphous 
bond-bending 3D network. 

While early measurements also identified a 
rigidity transition in these glasses near (r) = 
2.4, recent experiments have shown that there 
is not one but two transitions 0, 0, H, 0, El El 

mmmiiiQininins- startin s at iow 

mean coordination, the first transition seems 
to go from a floppy to a rigid but unstressed 
phase, and the second one on to a rigid and 
stressed phase. This new rigid but unstressed 
phase is known as the intermediate phase and 
is believed to be related to the self-organization 
of the system in order to minimize the stress in 
the network. 
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This interpretation is supported by a num- 
ber of theoretical works. In their pioneer- 
ing work, Thorpe and collaborators demon- 
strated that a network constructed by adding 
bonds with no stress allowed until it becomes 
inevitable would go through three phases: a 
floppy, a rigid-unstressed and a rigid-stressed 
phases 3J. The first transition, between the 
floppy and the rigid-unstressed phases, is the 
rigidity transition] the second one, between 
the rigid- unstressed and rigid-stressed phases, 
is the stress transition and happens immedi- 
ately when avoiding stress is no longer possible, 
The intermediate phase survives when the rigid 
but stress-free networks are fully equilibrated 
as was demonstrated by Barre and co-workers 
on Bcthe lattices |24| and Chubynsky et al. |30j 
on two-dimensional triangular lattices. This is 
not the only approach to self-organization, how- 
ever, and Micoulaut and collaborators [22, l23j | 
have shown that it is possible to recover an in- 
termediate phase in a stressed network if this 
stress is localized. 



III. METHODOLOGY 

The model we study here is the same as 
the one used in our previous paper |2jJ. I n 
our simulations, we use the pebble game al- 
gorithm described in the previous section; our 
computer code is based on the original program 
by D.J. Jacobs and M.F. Thorpe. 

We consider 2D triangular bond-diluted 
central-force networks. While we cannot make 
a direct comparison with experiment, previ- 
ous work has shown that the triangular lattice 
presents an intermediate phase similar to that 
of covalent glasses with angular constraints; our 
results should therefore be applicable, at least 
qualitatively, to experiments. 

In the original model of Thorpe and col- 
laborators |3j, bonds were added one by one 
and each checked for redundancy; redundant 
(stress-causing) bonds were rejected. Each new 
added bond was frozen in the network and 
was never moved nor removed. This proce- 
dure does not guarantee that stress-free net- 
works are equiprobable, making some networks 
more likely than others. To eliminate this bias, 
we introduced a bond-equilibration scheme, al- 
lowing the system to rearrange itself by moving 
bonds around. Each time a new bond is added, 
bonds are reshuffled throughout the lattice: a 
bond is picked at random, removed, and then a 
new bond is inserted in a random place choosing 
among those where it would not create stress. 
This bond-shuffling procedure is repeated un- 



til the system is equilibrated. We find that an 
equilibration of 10 iterations per added bond 
below (r}=3.5 and 100 iterations above (r)=3.5 
is enough for convergence (for more details, see 
Chubynsky et al. |30j) and it is the equilibra- 
tion scheme we use throughout this paper, un- 
less stated otherwise. In equilibrium, all stress- 
free networks with a given number of bonds (or 
mean coordination) are equiprobable. This cor- 
responds to the thermodynamic equilibrium at 
T — > for any model in which all stress-free 
networks have equal energy, but the energies of 
stressed networks are higher. 

As discussed in our previous paper |30j] . 
the intermediate phase in the model described 
above is associated with a non-trivial probabil- 
ity of finding a percolating network in this equi- 
librated model: this probability rises linearly 
from zero at the rigidity transition ((r) ~ 3.945) 
to one at the stress transition ((r) = 4.0), a 
result similar to that obtained by Barre et al. 
on the Bethe lattice |24|. Thus there are both 
percolating and non-percolating networks and 
these two classes need to be studied separately. 

All results presented in this paper (with the 
exception of the cluster size distributions given 
in Fig. |2l are obtained by running 200 indepen- 
dent simulations and obtaining the quantities of 
interest at different (r). When overall averages 
are presented, these are obtained by averaging 
over all these simulations. When, e.g., an aver- 
age for percolating networks is presented, then 
the averaging is done over only those of these 
networks that are percolating at the given (7'}. 
Obviously, close to the rigidity transition, very 
few of the 200 networks are percolating, and 
so the corresponding quantity will be an aver- 
age over a very small number of realizations and 
may contain a bigger error. When we report the 
results of insertion or removal of a single bond, 
only one attempt of insertion/removal per net- 
work is made, unless stated otherwise. 



IV. PROPERTIES OF RIGID 
CLUSTERS 

In this section, we look at properties of rigid 
clusters (both percolating and non-percolating) 
that exist in self-organized networks in the in- 
termediate phase. Some properties of rigid clus- 
ters that distinguish them from, e.g., clusters in 
usual connectivity percolation, need to be kept 
in mind. First, unlike in the connectivity case, 
a site can belong to more than one cluster (it 
then serves as a pivot joint between the clusters 
sharing this site). On the other hand, a bond in 
2D always belongs to just one cluster (not so in 
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3D, when it can serve as a hinge around which 
several clusters can rotate). For this reason, 
cluster decomposition in 2D is best expressed 
in terms of bonds and not sites. While this is 
not generally so, in the case of self-organized 
networks in the floppy or intermediate phase, 
the conversion between cluster sizes expressed 
in terms of bonds or sites is easy: since there 
is no stress, there is also no redundancy, and 
every rigid cluster of n sites contains exactly 
2n — 3 bonds. 

We start with properties of the percolating 
cluster. 
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A. Definition of a percolating cluster 

Normally, when the probability of percola- 
tion is either zero or one in the thermodynamic 
limit, the exact definition of percolation does 
not matter. In our case, since the probability of 
finding a percolating cluster increases linearly 
with average coordination in the intermediate 
phase, it is not clear that we can be so cavalier. 

For example, the probability that a cluster 
that percolates in only one direction (no matter 
which one) seems to be non-zero in the inter- 
mediate phase, although it is very low, around 
0.1 or less. This means that the results will dif- 
fer slightly if we define percolation using only 
one direction or both. To avoid confusion, wc 
choose here to call percolating those networks 
in which percolation occurs in both directions, 
and non-percolating those networks in which 
there is no percolation in either direction. Net- 
works with percolation in just one direction 
are ignored whenever we separate our results 
into those for percolating networks and those 
for non-percolating ones, but such networks are 
taken into account when this separation is not 
done. 



B. Size of the percolating cluster 

While we have already studied the probabil- 
ity that a percolating cluster occurs in the in- 
termediate phase, its size had not been charac- 
terized. Figure ^ shows the average number of 
bonds in the percolating cluster as a fraction of 
bonds actually present in the network. The av- 
eraging is done over all networks in which per- 
colation occurs. A remarkable feature is that 
even at the lowest (r) at which percolation is 
still (rarely) observed, the size of the percolat- 
ing cluster is well above zero. The smallest clus- 
ter size observed at the onset of the intermedi- 
ate phase is around 40%. This behavior is dif- 



FIG. 1: The fraction of bonds belonging to the per- 
colating rigid cluster among all bonds in the net- 
work, averaged over all percolating self-organized 
networks, for different network sizes indicated in 
the figure. All sizes here and in other figures in the 
paper are given in terms of sites. 



ferent from that for both connectivity and rigid- 
ity percolation on regular lattices in the random 
(non-self-organized) case, where the size of the 
percolating cluster, considered a good order pa- 
rameter, grows from zero at the transition, as 
expected in a second-order transition. This re- 
sult in the self-organized case is reminiscent of 
the first-order rigidity transition, such as that 
observed on Bethe lattices. 0, E^J 13^ • How- 
ever, our other results, as discussed below, do 
not support this analogy. 

Figure Q presents average sizes of perco- 
lating clusters. Given that the very exis- 
tence of the percolating cluster is uncertain in 
the intermediate phase (since only some net- 
works are percolating), it is reasonable to ask 
about the variation of the percolating cluster 
size. The quantity we look at is the stan- 
dard deviation, or width, of the fraction of 
bonds in the percolating cluster calculated as 
^/((F 2 ) - (F) 2 )n/(n - 1), where F is the frac- 
tion of bonds in the percolating cluster, (...) 
denotes the average over percolating networks 
and n is the number of percolating networks. 
In ordinary percolation, this width, of course, 
tends to zero as the network size grows; the 
percolating cluster size is a self-averaging quan- 
tity. Figure [3 shows that this is not so in our 
case. The width is above zero and is essentially 
size- independent. This is yet another difference 
from non-self-organized percolation (including 
that on Bethe lattices). Note at the same time 
that the width of the distribution of percolating 
cluster sizes is much smaller than the average 
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Mean coordination (r) 

FIG. 2: The standard deviation of the fraction of 
bonds in the percolating cluster for different net- 
work sizes, for self-organized networks. 

size. In other words, an overwhelming majority 
of networks have either a big percolating cluster 
or no percolating cluster at all — there are few 
(if any) "intermediate cases" with small perco- 
lating clusters. 



C. Sizes of non-percolating rigid clusters 

To further characterize the intermediate 
phase, it is useful to look at the distribution 
of rigid cluster sizes. In case of a second-order 
phase transition, the correlation radius is fi- 
nite away from the transition. As a result, at 
a certain cluster size there is a crossover from 
a power-law behavior to an exponential behav- 
ior. The correlation radius diverges as the tran- 
sition is approached from either side, and so 
the crossover moves towards bigger sizes as the 
transition is approached, and exactly at the 
transition, the power law persists indefinitely. 
Since the divergence of the correlation radius 
(or the crossover point) is governed by the same 
critical exponent on both sides of the transi- 
tion, it is expected that for two values of (r) 
at the same distance from the transition but on 
opposite sides, the crossover points will be the 
same by order of magnitude, if the transition 
is second-order. In case of a first-order transi- 
tion, the correlation radius does not diverge at 
the transition and there is always a crossover to 
the exponential behavior. 

Figure |21 shows the distribution of rigid clus- 
ter sizes at the coordinations of (r) = 3.92 (be- 
low the transition) and (r) = 3.97 (above the 
transition). The points are chosen at about the 
same distance from the rigidity transition, far 



enough from it to make sure that the correlation 
radius is not too big (if the transition is second- 
order) yet not too far to ensure that the points 
are still within the critical region. All results 
are obtained by running 200 independent simu- 
lations on networks of 40000 sites, each starting 
from an already equilibrated network and con- 
tinuing for 1000 additional equilibration steps. 
After each equilibration step, the distribution of 
non-percolating cluster sizes is obtained. While 
networks one equilibration step apart cannot be 
considered truly independent, inserting or re- 
moving a single bond often changes the rigid- 
ity of the network very significantly, justifying 
inclusion of the data obtained at every step. 
For (r) = 3.97, the distributions for percolating 
and non-percolating networks are plotted sepa- 
rately (at (r) = 3.92, obviously, there are very 
few percolating networks and we plot the re- 
sult for non-percolating networks only) . To de- 
crease the noise in the tail, all clusters of a given 
class (i.e., percolating and non-percolating) are 
binned using a logarithmic scale. 

Below the rigidity transition, in the floppy 
phase [(r) = 3.92; Fig. 01(a)], there is a clear 
crossover between the power-law and exponen- 
tial behaviors. The data are fitted using the 
product of a power law and an exponential: 

f x {x) = expt-x/xo). (3) 

We use the data between 10 2 and « 10 4 23 , 
i.e., dropping just a few data points at the tail, 
where noise and finite-size effects are signifi- 
cant, and omitting a region at the low end, as 
there are big deviations from the behavior de- 
scribed by Eq.|21 probably due to discretization 
effects from the lattice. The best fit is obtained 
with Ci = 3000, a x = 1.94 and x = 3900. The 
dashed line is a power law with the same C\ 
and oti, but without the exponential factor. 

Above the transition [(r) = 3.97; Fig. 02(b)], 
however, the power-law behavior persists at 
least for non-percolating networks, with no hint 
of the exponential tail, even for the largest sizes 
for which the data are available (around 30000, 
well above the crossover observed for (r) = 3.92 
around xo — 3900). For percolating nets, there 
is some deviation from the power law near the 
end, but it is likely due to finite-size effects 
(there is a percolating cluster taking up most 
of the network, so only relatively small non- 
percolating clusters are possible). To fit the 
data, we use pure power-law functions: 

f2{n, P } = C 2{ n, P }X- aHn * } , (4) 

where subscripts {n,p} refer to non-percolating 
and percolating networks, respectively. We use 



6 



10 4 


I 


i 1 i 


10 2 




■ 


10° 




- 


10" 2 


- 




io- 4 


- 




10' 6 


- \ d ) 

1 1 1 


1 — i — | — i — i-A 


10 




• non-perc. _ 


10 2 


- 


» perc. 


10° 






10' 2 






10' 4 






10' 6 


- (b) 

i 


1 , 1 \ 



1 10 100 1000 10000 



Cluster size 



FIG. 3: The distribution of sizes of non-percolating 
rigid clusters in non-percolating self-organized net- 
works at (r) = 3.92 [panel (a)] and in both perco- 
lating and non-percolating self-organized networks 
at (r) = 3.97 [panel (b)]. The details of the simu- 
lation and the fits (lines) are given in the text. 



data above 10 2 in the non-percolating case and 
from 10 18 to 10 28 in the percolating case. The 
values of the parameters providing the best fits 
are: C 2n = 6800, a 2n = 2.12, C 2p = 9100, 
a 2p = 2.74. While the difference between ot\ 
and a 2n is probably due to finite size and sam- 
pling effects, the difference between a 2n and 
a 2p suggests that these two quantities are dif- 
ferent, an unusual behavior. 

Our results indicate that the power-law dis- 
tribution of cluster sizes is observed in the 
whole intermediate phase, rather than at a sin- 
gle point (as would be the case without self- 
organization). In effect, the self-organization, 
which minimizes the stress in the network, 
maintains the system in a critical state through- 
out the intermediate phase. Similarly, unlike in 
the case of the critical point of random percola- 
tion, the cluster size critical exponent appears 
to be different in the non-percolating and per- 
colating cases. 

Based on the results of this section, we see 



that the rigidity phase transition in our model 
of self-organized percolation is very different 
from both first- and second-order phase tran- 
sitions. The minimization of the stress appears 
to act similarly to an external driving force, 
leading to a self-organized critical phase in a 
thcrmodynamically equilibrated system. 



V. RESPONSE TO A LOCAL 
PERTURBATION 

Given that both percolating and non- 
percolating networks coexist in the intermedi- 
ate phase, it is interesting to investigate the re- 
lation between the two classes. As we show 
in this section, a local perturbation involving 
the addition or removal of even a single bond is 
enough to affect rigidity of huge regions of the 
network and often converts a non-percolating 
network into a percolating one and vice versa. 

This behavior is not observed in regular ran- 
dom networks: since the probability of perco- 
lation in these systems is always either zero 
or one in the thermodynamic limit, a single 
bond can only change the percolation property 
right at the transition. Since a percolating net- 
work at that point is fractal and involves only 
an infinitesimal fraction of bonds, if an infinite 
cluster is created or destroyed, this can only 
involve an infinitesimal fraction of bonds and 
sites; away from the transition, the size of the 
affected region is always finite and then the 
fraction is obviously infinitesimally small. 

A few general comments about the conse- 
quences of addition or removal of a single bond 
are in order. First of all, the removal of a single 
bond can only break up the cluster to which this 
bond belongs; other clusters are not affected. 
This is because in 2D, all rigid clusters are al- 
ways rigid by themselves, i.e., they remain rigid 
when taken in isolation from the rest of the net- 
work. Conversely, the addition of a bond can 
merge several rigid clusters into one, but will 
not affect the clusters outside the resulting clus- 
ter. Note that even though the self-organized 
networks are by definition stress-free, a newly 
inserted bond can be redundant and introduce 
stress (always confined to the cluster in which 
it is inserted); in this case the configuration of 
rigid clusters is not affected, but the created 
stressed region may still be macroscopic and 
percolate. While we will briefly consider this 
situation at the end of this section, we mostly 
concentrate on the case when the inserted bond 
is non-redundant and thus affects the rigidity 
of the network but does not create stress. For 
brevity, we will call the places where insertion of 
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FIG. 4: The probability that a bond added to a 
non-percolating self-organized network at a random 
allowed place makes this network percolating. 



FIG. 5: The fraction of non-percolating self- 
organized networks that can become percolating 
upon a single bond addition. 



a bond does not create stress, as well as a bond 
inserted in such a place, allowed. Note that if 
we add an allowed bond creating a certain rigid 
cluster and then remove this bond and insert 
another allowed bond both ends of which are 
in the region coinciding with the cluster cre- 
ated by the first bond, exactly the same cluster 
will be created: indeed, the count of constraints 
within the region will be the same (2n — 3 for 
a region of n sites) , regardless of where in this 
region the bond is inserted. 



A. Conversion between percolating and 
non-percolating networks 

We first look at the conversion between non- 
percolating and percolating networks as a result 
of bond addition and removal. There are sev- 
eral related questions here. 



1. Addition of a bond in an allowed position 

In Fig. 0] we plot the probability that in- 
serting a bond at a randomly chosen allowed 
place makes a non-percolating network perco- 
lating. In the floppy phase (below (r) « 3.945), 
this probability tends to zero in the thermody- 
namics limit; however, it is non-zero everywhere 
in the intermediate phase. This probability is 
close to zero just above the rigidity transition 
and approaches one close to the upper bound- 
ary of the intermediate phase at (r) = 4. 

We can also ask about the frequency of net- 
works in which conversion from non-percolating 
to percolating upon a single bond addition is 



possible. Note that because of the above- 
mentioned property that insertion of a bond 
into a region made rigid by another bond rigid- 
ifies exactly the same region, we do not need to 
try each and every allowed bond, which speeds 
up the simulation enormously. Fig. |3] plots the 
fraction of non-percolating networks that can 
become percolating with a bond placed judi- 
ciously in an allowed position. Note that this 
is different from the previous figure: instead of 
choosing a place to insert a bond at random, we 
now make the best effort to cause percolation 
with a single bond addition, if at all possible. 
As is expected, this quantity is higher than that 
plotted in Fig. 0] but the difference is small; 
over most of the intermediate phase, whenever 
there are any allowed positions where bond in- 
sertion creates a percolating cluster, most al- 
lowed positions will do. 

It is also interesting to note that the quantity 
in Fig.[5]is very close to linear in the intermedi- 
ate phase, and is probably exactly linear, just 
like the probability of percolation without any 
bond insertions. In the Appendix, we explain 
why these two quantities are equal. 

Finally, in Fig. the average size of the per- 
colating cluster arising after bond insertion is 
compared to the average size of the percolat- 
ing cluster in those cases when it exists even 
without bond insertion (i.e., the quantity in 
Fig. These values are identical: the size 
of the percolating cluster emerging after insert- 
ing a bond in a non-percolating network is the 
same as in originally percolating networks. This 
suggests that a percolating network that arises 
after bond insertion is a typical percolating net- 
work, just like those networks that percolate 



s 



S 

o 



o 



O 
CO 



0.8 



0.6 



0.4 




without insertion 
after insertion 



3.94 3.96 3.98 

Mean coordination (r) 



FIG. 6: The comparison between the average size 
of the percolating cluster in originally percolating 
self-organized networks ("without insertion") and 
in originally non-percolating networks that become 
percolating after bond insertion ("after insertion"). 
The latter quantity is calculated as the average over 
all networks that can become percolating, by using 
"judicious placement" of a bond, as described in 
the text. 



"without insertion. In the next subsection, wc 
present more evidence in favor of this. Also see 
a more detailed discussion of this and a possible 
caveat in the Appendix to this paper. 

Given the significant probability of conver- 
sion of non-percolating networks to percolating 
ones (see Fig. |S} , we can hypothesize that all 
non-percolating networks can become percolat- 
ing after a finite number of bond insertions; 
moreover, the average size of the percolating 
cluster after the minimal number of insertions 
needed to create it is again the same as the 
average size of the percolating cluster in net- 
works that are percolating without bond inser- 
tions. An even stronger hypothesis is that af- 
ter each new bond addition, the fraction of so 
far non-percolating networks becoming perco- 
lating is still the same linear function of (r) as 
in Fig. [5] These hypotheses need to be tested 
in the future. 



2. Removing a bond from the network 

Likewise, we can consider the probability 
that removing a randomly chosen bond from 
a percolating network breaks the percolating 
cluster. This quantity is shown in Fig. Just 
like for bond removal, this probability is non- 
zero everywhere in the intermediate phase. In- 
terestingly, even at the upper boundary ((r) — > 
4), where almost all networks are percolating, 



it is still very easy to break percolation and 
thus create a non-percolating network. In fact, 
the probability is the highest in this limit. Of 
course, this can be attributed to the fact that 
the percolating cluster is itself the biggest at 
this point (taking up the whole network), so 
there is a greater chance than elsewhere to se- 
lect a bond that belongs to it (which is, of 
course, a necessary condition of its destruction) . 
This effect can be factored out by dividing the 
quantity in Fig. by the average fraction of 
bonds in the percolating cluster. This will 
then give the probability of the destruction of 
the percolating cluster, given that the removed 
bond belongs to this cluster. This is plotted in 
Fig. [S] It is seen that across the whole interme- 
diate phase, this quantity is nearly constant at 
around 70%. In other words, everywhere in the 
intermediate phase, the percolating cluster on 
average contains about 70% of bonds such that 
removal of any one of them will destroy percola- 
tion. This is true even close to (r) = 4. This is, 
of course, just the average; one could ask if for 
some networks this quantity is zero (similar to 
how in the case of bond insertion, not all non- 
percolating networks can be made percolating 
by a single bond addition, as shown in Fig. [5}. 
To check if this is the case, we remove (and 
then reinsert) up to 10 bonds one by one (all 
chosen within the percolating cluster) and see 
if the percolating cluster ever gets destroyed. 
This happened for 5398 out of 5429 90000-site 
networks, or about 99.4%, including 198 out of 
200 networks (99%) at (r) = 3.999 (the highest 
mean coordination in our simulations). Since 
these fractions are so close to 100%, we can hy- 
pothesize that, unlike in the case of bond inser- 
tion, in fact, all percolating networks have some 
bonds whose removal destroys percolation. 

To conclude this subsection, our results for 
the change in the percolation state of the net- 
work upon single bond addition or removal 
again confirm that the system remains critical 
in the whole intermediate phase. 



B. Change of rigidity: "rigid" and 
"floppy" bonds 

In the previous subsection, we have found 
that in many cases insertion or removal of a 
single bond can change the percolation status 
of the network and affect the rigidity of its sig- 
nificant part. Since we were dealing with perco- 
lating clusters only, we could not study the ef- 
fect of bond insertion or removal in those cases 
when the percolation status does not change. 
To do this, we introduce the concept of "rigid" 
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FIG. 7: The probability that removal of a randomly 
chosen bond from a percolating self-organized net- 
work destroys percolation. 
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FIG. 8: The probability that removal of a random 
bond destroys percolation (the quantity in Fig. [7f 
divided by the average fraction of bonds in the per- 
colating cluster (the quantity in Fig.0. This serves 
as an estimate of the probability that a bond chosen 
at random among those belonging to the percolat- 
ing cluster destroys percolation. 



and "floppy" bonds. 

Below the rigidity transition, only small rigid 
clusters are present in the network. Interest- 
ingly, a significant fraction of them consist of 
just a single bond. The number of bonds be- 
longing to these single-bond clusters decreases 
as the rigidity transition is approached and 
eventually crossed, although they are still en- 
countered well above the transition in floppy 
pockets of the network. Since these bonds are 
associated with the floppy phase and floppy re- 
gions of the network in the rigid phase, we call 
such bonds floppy. All other bonds (i.e., those 
belonging to clusters consisting of more than 



FIG. 9: The average fraction of rigid bonds in the 
network in the random and self-organized cases (in 
the latter case, also separately for percolating and 
non-percolating networks). All simulations are for 
networks of 90000 sites, the overall averages are 
over 200 networks, the averages restricted to perco- 
lating and non-percolating networks are over those 
of the 200 self-organized networks that are respec- 
tively percolating and non-percolating. 



one bond, or, in other words, rigid with respect 
to some other bonds) are called rigid. The frac- 
tion of rigid bonds as a function of (r) is plot- 
ted in Fig.|U|for both random and self-organized 
networks. As expected, it grows in both cases 
with increasing (r), but in the self-organized 
case it reaches 1 at (r) = 4 (when the network 
becomes fully rigid), which in the random case 
does not happen until the full coordination at 
(r) = 6. For the self-organized case, the av- 
erages over just percolating networks and over 
just non-percolating networks are also shown 
in Fig. |S| Interestingly, the average for non- 
percolating networks remains nearly constant 
over the whole intermediate phase at about 
75%, even as (r) — > 4. In this limit, the num- 
ber of floppy modes per site tends to zero in all 
networks; yet, as these results show, in those 
few networks that still do not percolate, many 
bonds (about 25%) are still in single-bond clus- 
ters. 

Just as in the previous section we considered 
the width of the distribution of percolating clus- 
ter sizes, it is interesting to look at the width 
of the distribution of fractions of rigid bonds. 
Of course, since the averages are very differ- 
ent for percolating and non-percolating net- 
works, it makes sense to separate these two 
classes. The results are shown in Fig. Note 
that in the intermediate phase, the width is 
size-independent and so likely remains finite 
in the thermodynamic limit, for both perco- 
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FIG. 10: The standard deviation of the distribution 
of fractions of rigid bonds in percolating [panel(a)] 
and non-percolating [panel (b)] self-organized net- 
works for different network sizes. 



lating and non-percolating networks, just as 
we have seen for the percolating cluster size 
(see Fig. 0). But in the floppy phase (where, 
of course, only non-percolating networks are 
present), the width clearly decreases fast with 
size (it is roughly inversely proportional to the 
square root of the network size) . Thus the num- 
ber of rigid bonds is a self-averaging quantity 
in the floppy phase but not in the interme- 
diate phase, even when percolating and non- 
percolating networks are considered separately. 
At the same time, the widths for both percolat- 
ing and non-percolating cases are much smaller 
that the difference between these two cases, so 
based on the count of rigid bonds, these two 
classes are clearly distinct. 

We now look at the influence of insertion or 
removal of a bond on the number of rigid bonds 
and their spatial distribution. 

In Fig. ^] we show the average change in 
the number of rigid bonds upon insertion of a 
single bond, in both the random and the self- 
organized cases. Note that in the random case, 
this change is very small on average, around 
10 bonds or less. On the other hand, in the 
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FIG. 11: The average change in the number of rigid 
bonds upon addition of a single bond, for random 
and self-organized networks of 90000 sites. 



self-organized case this quantity diverges very 
fast when the rigidity transition is approached. 
This is not surprising: we have seen that infinite 
percolating clusters can easily form and break 
in this case. For this reason, it makes sense to 
look at the fraction of bonds that undergo the 
change. 

In Fig. El we show the fraction of bonds 
converted from floppy to rigid when a bond is 
added to the network at a randomly-chosen al- 
lowed position, among all floppy bonds in the 
network. In panel (a), we plot the overall av- 
erage, as well as partial averages restricted to 
those cases when the network goes from non- 
percolating to percolating, remains percolating 
and remains non-percolating (the curve for the 
latter case is barely above the baseline) . In pan- 
els (b), (c), and (d), we show these partial aver- 
ages for several different sizes. Based on what 
we have already learned for networks switching 
from non-percolating to percolating, it is quite 
natural that in this case a significant fraction 
of floppy bonds become rigid; this fraction ap- 
proaches 1 when (r) — > 4, which is again ex- 
pected, since in this limit the percolating clus- 
ter takes up the whole network (see Fig. [fj). 
The other two cases are more interesting. In 
the floppy phase, only the non-percolating — > 
non-percolating situation is possible, and as we 
see from panel (d), in this region of the phase di- 
agram the average fraction of bonds converted 
from floppy to rigid decreases fast with size. In 
the inset, we plot the average number (rather 
than fraction) of converted bonds, and we see 
that this number is size-independent. So in the 
floppy phase, a finite number of bonds gets con- 
verted. In the intermediate phase, the situation 
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is different. Looking first at the percolating 
— > percolating situation [panel (c)], the frac- 
tion of converted bonds depends only weakly 
on the network size, and it is possible that this 
quantity goes to a constant in the thermody- 
namic limit. Even if, in fact, this fraction de- 
creases to zero as the size N goes to infinity, it 
is clear that the decrease is much slower than 
oc 1/N, and thus the mean number of converted 
bonds diverges when N — > oo. In the non- 
percolating — > non-percolating case, the quality 
of the data is lower (there is very much noise 
due to a very large variation in the number of 
converted bonds), but still it is clear that the 
decrease with N (if present) is certainly much 
slower than oc 1/N. 

The fact that for the percolating — > perco- 
lating and non-percolating — ► non-percolating 
cases the average number of converted bonds di- 
verges in the thermodynamic limit in the inter- 
mediate phase means that even when the perco- 
lation status does not change, the region of the 
network whose rigidity is affected are macro- 
scopic at least in some cases. This is again 
consistent with the criticality of the interme- 
diate phase. The distribution of the sizes of af- 
fected regions (or numbers of converted bonds) 
is likely power-law, which needs to be tested in 
the future. 

As an illustration of effects of bond insertions 
we show two examples in Fig. 1131 The up- 
per panel shows an example for the case when 
the network switches from non-percolating to 
percolating after a bond is added. The lower 
panel shows an event where the network re- 
mains non-percolating, but large-scale rigidifi- 
cation still occurs without percolation. In both 
cases, the added bond is red and pointed with 
an arrow, thick green bonds are those that are 
originally rigid (and, of course, remain rigid af- 
ter bond addition), thin blue bonds are orig- 
inally floppy and remain floppy, finally, thick 
black bonds are of most interest: these are the 
ones that switch from floppy to rigid. In the 
first (non-percolating — * percolating) case, con- 
verted bonds are spread throughout the net- 
work. Many bigger rigid clusters separated by 
floppy "interfaces" merge into one percolating 
rigid cluster; in essence, the figure illustrates 
the rigidification of these interfaces. In the sec- 
ond (non-percolating — > non-percolating) ex- 
ample, the affected region is still large, but non- 
percolating. This example is larger than aver- 
age, as can be deduced from Fig. El panel (d), 
but it is not a rare event. 

We find overall similar behavior upon re- 
moval of a bond. The results for the fraction of 
rigid bonds converting to floppy are in Fig. EI 



Similarly to Fig. El in panel (a) we have the 
overall average and partial averages for perco- 
lating — > non-percolating, percolating — > perco- 
lating and non-percolating — ► non-percolating 
cases, for a single network size (90000 sites). 
In panels (b), (c), and (d), we have the same 
partial averages, but for three different sizes. 
Again, the conclusions are similar to the case 
of bond insertion: in the percolating — * non- 
percolating case, when the percolation status 
of the network changes [panel (b)], the fraction 
of bonds that switch from rigid to floppy is ex- 
pectedly high. It is much lower in the other 
two cases, when the percolation status does 
not change, but still, just as for bond addition, 
while in the floppy phase the fraction of con- 
verting bonds falls rapidly with increasing size 
and, as the inset of panel (d) shows, the aver- 
age number of converting bonds remains con- 
stant, in the intermediate phase, again, the de- 
pendence of the fraction of converting bonds on 
the size is very slow and the number of convert- 
ing bonds diverges in the thermodynamic limit 
- thus again, macroscopic regions of the net- 
work can be involved. 

We finish this subsection with an observa- 
tion: there exists a symmetry in bond conver- 
sions upon addition and upon removal. Namely, 
the average number of bonds converting upon 
bond addition in the case when the network 
transforms from non-percolating to percolating 
is the same as the average number of bonds con- 
verting upon bond removal when the network 
transforms from percolating to non-percolating. 
This is illustrated in Fig. El in this figure, 
we plot the ratios of the above-mentioned num- 
bers and the total number of bonds in the net- 
work. We find that these two quantities coin- 
cide. Note that these quantities are different 
from those plotted in Figs. 1121 and!14l in these 
figures, the numbers of converted bonds were 
divided by the number of floppy bonds and the 
number of rigid bonds, respectively, and not by 
the total number of bonds. The equality is easy 
to understand, if we recall that we have already 
seen some evidence (see Fig. |HJ| that networks 
that become percolating after bond insertion 
are, in fact, typical percolating networks, just 
like those that are originally percolating. By 
extension, we can assume that networks that 
become non-percolating after bond removal are 
also typical non-percolating networks. If so, 
then the average change in the number of rigid 
bonds should in both cases be the same as the 
difference in the average number of rigid bonds 
between percolating and non-percolating net- 
works. Indeed, in Fig. El we also plot the dif- 
ference between the average fractions of rigid 
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FIG. 12: The average fraction of floppy bonds converting to rigid upon addition of a single bond to a self- 
organized network. Panel (a) shows the overall average, as well as partial averages for the cases when the 
network converts from non-percolating to percolating, remains non-percolating and remains percolating, 
for networks of 90000 sites. Panels (b), (c), and (d) show the partial averages in the non-percolating 
— > percolating, percolating — » percolating, and non-percolating — > non-percolating cases, respectively, for 
three different network sizes. The inset in panel (d) shows the average number of converted bonds in the 
non-percolating — > non-percolating case. 



bonds in the non-percolating and percolating 
cases (i.e., between the dashed and the dot- 
dashed lines in Fig. |§J) ; it is seen that this quan- 
tity coincides with the other two. 



C. Stress propagation 

In the previous subsections, we have looked 
at cases when an "allowed" bond insertion is 
done; in other words, the bond is inserted in 
one of those places where it does not create 
stress. We now look at the opposite situation, 
i.e., we analyze the results of inserting a bond in 
one of "disallowed" places. In this case, the in- 
serted bond is redundant, which means that its 
insertion does not change the configuration of 



rigid clusters, but some bonds (including the in- 
serted one) become stressed. We are interested 
in the emerging stressed region and in particu- 
lar, whether it percolates or not. 

Similarly to the case of allowed bond in- 
sertion, we first look at the probability that 
the stressed region percolates. Note that the 
stressed region emerging upon bond insertion 
is always restricted to the rigid cluster contain- 
ing the new bond. For this reason, the stressed 
region can only percolate if the original network 
is percolating (i.e., contains a percolating rigid 
cluster) , and we only need to look at percolating 
networks. In Fig. 1161 we show the probability 
that in a network in which rigidity percolates, 
a percolating stressed region emerges after in- 
sertion of a bond at a randomly chosen disal- 
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the rigidity percolation status of the network. 




Again, we see that the size of the stressed re- 
gion as a finite fraction of the total number of 
bonds is only weakly size-independent and thus 
the affected region can be macroscopic in the 
intermediate phase even when there is no perco- 
lation. In the floppy phase, on the other hand, 
the fraction of stressed bonds decays fast with 
increasing network size and as the inset (where 
the number of conversions, rather than the frac- 
tion, is plotted) shows, the average number of 
conversions is size-independent and the affected 
region remains finite. 




FIG. 13: Two examples of rigidification of self- 
organized networks upon addition of a bond. In 
the top panel, the network switches from non- 
percolating to percolating. In the bottom panel, it 
remains non-percolating. In both cases, the added 
bond is red and pointed with an arrow, thick green 
bonds are those that are originally rigid, thin blue 
bonds remain floppy, thick black bonds switch from 
floppy to rigid. Both networks contain 10000 sites. 

lowed place (not necessarily within the perco- 
lating cluster). We note that apparently, this 
probability remains finite (does not go to zero) 
in the thermodynamic limit. 

We also look at the sizes of stressed regions 
arising upon disallowed bond insertion. The re- 
sults are in Fig. El Panel (a) shows the overall 
average fraction of stressed bonds and the par- 
tial averages for cases when stress does and does 
not percolate, for networks of 90000 sites; panel 
(b) gives the partial average for the case when 
stress percolates; and panel (c) gives the partial 
average for the case when stress does not per- 
colate. We do not make distinctions based on 



VI. EXPERIMENTAL EVIDENCE FOR 
THE SELF-ORGANIZED CRITICAL 
BEHAVIOR 



An atypical response to an external perturba- 
tion seen in the intermediate phase has been ob- 
served experimentally [Tsl ] . Raman pressure ex- 
periments show that applying external pressure 
to a network blueshifts the frequency of corner- 
sharing (CS) tetrahedral units only once a cer- 
tain pressure threshold (P c ) value is reached. 
P c is found to be zero within the intermedi- 
ate phase and non-zero outside the intermediate 
phase. These experimental results, published in 
Ref. 0| , are reproduced in Fig.HTfl 

Outside the intermediate phase, the presence 
of a non-zero threshold for the shift in the 
frequency suggests that inhomogeneities can- 
cel the effect, leaving only a broadening of the 
peak [3. On the other hand, we have seen that 
in the intermediate phase, a finite fraction of 
the network is either in a percolating cluster or 
in a region that becomes percolating after a few 
bond insertions. This region is also stress-free. 
This can only be achieved by a precise balanc- 
ing between constraints and degrees of freedom 
not just on average, but on all length scales, 
in other words, by the network being homoge- 
neous because of self-organized criticality. We 
note that the way rigidity and stress propagate 
through a macroscopic fraction of the network 
in our simulations is analogous to how stress 
propagates uniformly when external pressure is 
applied in experiments. While the relation be- 
tween P c and inhomogencity still needs to be 
fully established, self-organized criticality in the 
intermediate phase is consistent with the obser- 
vations. 
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FIG. 14: The average fraction of rigid bonds converting to floppy upon removal of a single bond from a self- 
organized network. Panel (a) shows the overall average, as well as partial averages for the cases when the 
network converts from percolating to non-percolating, remains non-percolating and remains percolating, 
for networks of 90000 sites. Panels (b), (c), and (d) show the partial averages in the percolating — > 
non-percolating, percolating — > percolating, and non-percolating — > non-percolating cases, respectively, for 
three different network sizes. The inset in panel (d) shows the average number of converted bonds in the 
non-percolating — > non-percolating case. 



VII. CONCLUSION 

We have studied the structural and response 
properties of the intermediate phase in the 
phase diagram of rigidity percolation using a 
model of self-organization on a 2D triangular 
network. 

We had shown previously that the proba- 
bility of rigidity percolation in the intermedi- 
ate phase increases linearly from zero to one 
as a function of the mean coordination (r). 
At any (r), there are both percolating and 
non-percolating networks in the ensemble. In 
this paper, we have looked at the properties 
of both percolating and non-percolating clus- 
ters, the latter separately for percolating and 
non-percolating networks. It turns out that 



at the point at which the percolating networks 
first emerge (the rigidity transition), the perco- 
lating cluster takes up about 40% of the net- 
work, unlike the case of the usual second-order 
phase transition, where the emerging cluster is 
fractal and thus the fraction of bonds belong- 
ing to it is zero in the thermodynamic limit. 
The size of the percolating cluster and some 
other characteristics are not self-averaging, but 
the distributions are rather narrow, if perco- 
lating and non-percolating networks are con- 
sidered separately. The distribution of sizes of 
non-percolating clusters is exponential at large 
sizes in the floppy phase, but power-law for ar- 
bitrarily big cluster sizes in the whole inter- 
mediate phase and not just at the transition. 
The power-law exponents are different for non- 



15 



CO 
"D 

c 
o 



o 



0.3 



0.2 



5 0.1 




removal 
addition 
difference 



3.94 3.96 3.98 

Mean coordination (r) 



FIG. 15: The average fractions of bonds: convert- 
ing from rigid to floppy when a bond is removed 
and the network switches from percolating to non- 
percolating (green); converting from floppy to rigid 
when a bond is inserted and the network switches 
from non-percolating to percolating (red). These 
are compared to the difference between the average 
fractions of rigid bonds in the percolating and non- 
percolating cases. All quantities are for networks 
of 90000 sites. 
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FIG. 16: The probability that a percolating 
stressed region forms when a bond is inserted at 
a random "disallowed" place into a network with a 
percolating rigid cluster. 



as that for the initially percolating networks. 
In a sense, the percolating cluster is "hidden" 
in the non-percolating network and is revealed 
upon addition of a single bond. It appears, 
further, that all non-percolating networks can 
be converted into percolating ones with a finite 
number of bond additions, with the same size of 
the resulting cluster as in networks that perco- 
late without bond addition. This implies that 
in the thermodynamic limit, there is no differ- 
ence between percolating and non-percolating 
networks. The intermediate phase then can be 
thought of as being the region of the phase di- 
agram where all networks possess a percolat- 
ing region that is nearly isostatic. That is, 
this region can technically be floppy or stressed, 
but the number of floppy modes or redundant 
constraints causing stress is negligible in the 
thermodynamic limit. The lower boundary of 
the intermediate phase is then, strictly speak- 
ing, not the rigidity percolation transition, but 
rather the "near-isostaticity" percolation tran- 
sition. 

These unusual properties indicate that the 
intermediate phase is a self-organised critical 
phase, with the system staying at the rigid- 
ity percolation threshold for the whole range 
of mean coordinations. The existence of such a 
self-organized nearly isostatic critical phase ex- 
plains recent pressure experiments |l8j . Since 
rigidity percolation occurs on a finite fraction 
of the network under any microscropic pertur- 
bation, the strain associated with the external 
pressure is immediately transferred to a macro- 
scopic fraction of the sample, leading to a shift 
in the Raman spectrum. 

We can certainly expect other surprises asso- 
ciated with this intermediate phase in the rigid- 
ity phase diagram. But for that, we will likely 
need to create more realistic models of the ex- 
perimental systems. 



percolating and percolating networks. 

We have also looked at the changes in the 
rigidity of the network due to a microscopic per- 
turbation in the form of insertion or removal 
of a single bond. It turns out that one bond 
is often enough to convert a non-percolating 
network into a percolating one and vice versa; 
moreover, when a non-percolating network is 
turned into a percolating one, the average size 
of the resulting percolating cluster is the same 
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Appendix: Relation between the 
probability of percolation in the 
intermediate phase and the probability of 
inducing percolation by single bond 
insertion 



Numerically, the probability that there ex- 
ists a place in a non-percolating network such 
that insertion of a bond at that place causes 
percolation is the same linear function in the 
intermediate phase as the probability of rigid- 
ity percolation (see Fig. [SJ . Here we show that 
this equality is a necessary consequence of the 
latter quantity being finite (between zero and 
one) and continuous throughout the intermedi- 
ate phase. 

Given the probability of percolation among 
networks with B bonds, we can calculate this 
probability for networks with B + 1 bonds. 
The method is similar to that used to calcu- 
late the bond-configurational entropy of self- 
organized networks in our previous paper [30|. 
Suppose there are Nbc(B) networks (or bond 
configurations) with B bonds that are stress- 
free. If the probability of percolation is s(B), 
then N P (B) = s(B)Nh c (B) of these are perco- 
lating and the rest, N n (B) = [1 - s(B)]N hc (B), 
non-percolating. Also suppose there are on av- 
erage N ap (B) allowed places to insert a bond 
in percolating networks and N an (B) in non- 
percolating networks and the probability that 
a non-percolating network becomes percolat- 
ing upon a random bond insertion is tt(B). 
Since all percolating networks will remain per- 
colating upon bond insertion, then an aver- 
age percolating network with B bonds will 
produce N ap (B) different percolating networks. 
On the other hand, only a fraction tt(B) of 
non-percolating networks will become percolat- 
ing, so an average non-percolating network will 
produce Tr(B)N an (B) percolating networks and 
[1 — ir{B)]N an {B) non-percolating networks. If 
we simply multiply these numbers by the num- 
ber of networks with B bonds of each kind, this 
will not produce the correct count of percolat- 
ing and non-percolating networks with B + 1 
bonds, because each such network can be ob- 
tained in many different ways. Specifically, each 
network with B + 1 bonds can be produced by 
bond insertion from as many different networks 
with B bonds as there are the latter that can 
be obtained by bond removal from the former. 
This number is always B + 1, since removal of 
every bond will produce a distinct network and 
all of them are stress-free. So the count has to 
be divided by B + 1 , and we get for the number 



of percolating networks with B + 1 bonds, 



N p {B + l) = 
N p (B)N ap (B) 



N n (B)ir(B)N aD (B) 



B + l 



(5) 



and for the number of non-percolating networks 
with B+l bonds, 

, N n (B)\l - TT(B)]N an (B) /N 

N n (B + l)= " v n v n v ; . (6) 

B + l 

Then the probability of percolation for a net- 
work with B + l bonds is 



s(B + l) = 



N p (B + 1) 



N p (B + 1) + N n (B + 1) 
N p (B)N ap (B) + N n (B)ir(B)N an (B) 
N p (B)N ap {B) + N n (B)N an (B) 



(7) 



Note that since the probability of percolation 
changes continuously, the difference between 
s(B + 1) and s(B) is 0(1/N), where N is 
the network size. Neglecting terms that are 
0(1/N), we should equate s(B + 1) = s(B) 
and then, omitting the argument B for brevity 
and using N p / {N n + N p ) = s, 



sN ap + tt(1 - s)N a 
sN ap + (1 - s)N an 



N a 



(8) 



(9) 



ap 



Now, as we have seen, networks that become 
percolating after bond insertion are typical, 
generic percolating self-organized networks. So 
we can expect that on average, when a bond 
is inserted into a non-percolating network and 
it becomes percolating, the change in the num- 
ber of allowed bonds is Nan — N ap . But bonds 
that become disallowed are exactly those that 
are within the percolating cluster, and, accord- 
ing to our arguments at the beginning of Sec- 
tion these and only these bonds will cre- 
ate the percolating cluster and cause percola- 
tion (indeed, if insertion of any bond creates a 
percolating cluster, then insertion of any bond 
outside the region that turns into the percolat- 
ing cluster will not cause percolation because 
the cluster that will form cannot overlap with 
this region, other than over a single pre-existing 
cluster, and thus cannot percolate). So, if for a 
given network making it percolating after one 
bond insertion is at all possible, then the prob- 
ability that percolation will occur after ran- 
dom insertion is equal to the fraction of al- 
lowed bonds within the percolating cluster, or 
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(Nan — N ap )/N an , and then the right-hand side 
of Eq. {HI is exactly the probability that perco- 
lation after one bond insertion is possible, and 
the proof is complete. 

There are some caveats in this proof, related 
to the fact that, as we have seen, many quan- 
tities in the intermediate phase are not self- 
averaging. Because of this, for instance, if frac- 
tions of allowed bonds are not strictly deter- 
ministic, there can be correlations between the 
number of allowed bonds in a network and this 
number in the percolating network created af- 
ter bond insertion. For instance, the latter may 
tend to be above average when the former is 
below average. Then the average change in the 
number of allowed bonds may differ somewhat 
from N &n — iV ap , which is the difference of av- 
erages. In fact, the very statement that net- 
works obtained by bond insertion can be con- 
sidered typical needs to be treated with caution. 



Indeed, if all quantities associated with self- 
organized networks were self-averaging, then, in 
order to change the values of these quantities by 
inserting a bond, one would need to upset the 
balance so badly as to create a finite chance of 
seeing networks with properties that at equi- 
librium have probability zero. Such bias after 
inserting just a single bond cannot happen, ex- 
cept in very special cases. But the situation 
changes when there is no self-averaging: if the 
distribution of values of a certain quantity has 
a finite width, then any bias can change the av- 
erage. For this reason, it is possible that both 
the equality proved above and the statement 
that any quantities obtained for networks that 
become percolating after inserting a bond are 
the same on average as for originally percolat- 
ing networks (as demonstrated in Figs. El and 
I15fl are, in fact, only approximate. In any case, 
deviations (if any) are very small. 
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FIG. 17: The average fraction of stressed bonds 
in the self-organized network after a single "disal- 
lowed" bond is inserted. Panel (a) shows the over- 
all average, as well as the partial averages for cases 
when stress does and does not percolate, for net- 
works of 90000 sites. Panel (b) shows the partial 
average for the case when the network percolates, 
for three different network sizes. Panel (c) shows 
the partial average when the network does not per- 
colate, again, for three different network sizes; in 
the inset, the average number of stressed bonds is 
shown for the same case. 
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FIG. 18: Variations in the frequency of CS tetra- 
hedral units as a function of pressure for different 
Ge^Sei-^ glasses. Triangles are results taken from 
the work of Murase and Fukunaga |33| and filled 
circles are results of the work of Wang et al. |l8j. 
The figure is taken from Ref.fTj| 



